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ABSTRACT 

Using 3D hydrodynamic calculations we simulate formation of molecular clouds in 
the Galaxy. The simulations take into account molecular hydrogen chemical kinet- 
ics, cooling and heating processes. Comprehensive gravitational potential accounts for 
contributions from the stellar bulge, two and four armed spiral structure, stellar disk, 
dark halo and takes into account self-gravitation of the gaseous component. Gas clouds 
in our model form in the spiral arms due to shear and wiggle instabilities and turn into 
molecular clouds after t ^ 100 Myr. At the times t ~ 100 — 300 Myr the clouds form 
hierarchical structures and agglomerations with the sizes of 100 pc and greater. We 
analyze physical properties of the simulated clouds and find that synthetic statistical 
distributions like mass spectrum, "mass-size" relation and velocity dispersion are close 
to those observed in the Galaxy. The synthetic / — v (galactic longitude - radial veloc- 
ity) diagram of the simulated molecular gas distribution resembles observed one and 
displays a structure with appearance similar to Molecular Ring of the Galaxy. Exis- 
tence of this structure in our modelling can be explained by superposition of emission 
from the galactic bar and the spiral arms at ^3-4 kpc. 
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1 INTRODUCTION 

Molecular clouds represent one of the major con- 
stituents of our galaxy (Dame et al. 20011. Their masses 



fc Roberts||T98T 



Kwan & Valdes 



Tomisaka 1984; |Roberts fc Stewart|[T987l 



Solomon et al. ||1987| ). | Song | l999 l i Dobbs 2008) 



1987 \ or interaction of the turbulent flows 
(e.g. |Ballcstcr os^-Paredes et al. |2007 |. In reality both types 
of processes are expected to play significant role (Zhang & 



Molecular clouds play very important role in determination 
of the structure and evolution of the Galaxy. This can be 
illustrated by the fact that the stars and their clusters orig- 
inate from molecular clouds (Shu et al. 1987 McKee & 



Ostriker||2007| |Lada fc La da 2003). The dynamical time of 



disintegration for molecular clouds in the Milky Way is much 



shorter than the rotation period of the galaxy (e.g. Blitz & 
Shu|[l9 80 ) . So, explanation of the observational fact of ex- 
istence of the rich population of molecular clouds requires 
involvement of efficient mechanisms leading to formation of 
molecular clouds in the Galaxy. Special interest is paid to 
formation of the g iant molecular clou ds (GMC) which have 
masses Si 1O 6 M0 ( Dame et al. |l986 1 . Two main theories of 
the GMC formation were proposed. The first one explains 
formation of the GMCs by large scale magnetic and/or gravi- 



tational instabilities I Elmegreen 1979 Balbus & Cowie 1985 
Elmegreen||1994| |Chou et al. ||2000| |Kim fc Ostriker||2006 1 



In the second one the GMCs are built up by coalescence of 



the smaller molecular clouds ( Field & Saslaw 1965 Levinson 



Large scale perturbations in a gas of the Galactic disk 
can be generated not only by magnetic fields or self-gravity, 



but also by the galactic spiral shock wave (Roberts 1969 



Nelson & Matsuda 19771. A simple analysis of the global 
galactic shock wave shows that such shocks are unstable 
(|Mishurov &i Suchkov 19751. Obviously, large scale shear 



flows in the Galactic shocks lead to the development of both 
wiggle and Kelvin-Helmholz instabilities (Wada, 2001 Wada 



& Koda 2004 1. Such hydrodynamic instabilities are expected 



to give a birth to spurs, gaseous fragments and turbulent 



flows (Wada|1994 


Wada et al. |2000| 


Dobbs & Bonnel 2006 


Shetty & Ostriker 


20061. The gravitational and thermal (in 



general, thermo-chemical) instabilities are expected to de- 
velop further fragmentation, which leads to formation of the 
population of molecular clouds. 

The properties of molecular clouds (mass, size, density, 
temperature etc.) are intensively studied both observation- 
ally ( |Solomon et al. |1987||FaI garone et al. |1992||Lada et al. 
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2008; Hcyer et al. |2009 Rathborne et al. |2009; Kauffmann| 


et al. 


2010 Roman-Zuniga et al. 


2010 

1 TIT 


1 and numerically 



Tan|2009| |Glover fc Mac Low|2011| >. Efforts in order to un- 



derstand dependencies between physical parameters of the 
clouds lead to establishment of several empirical relation- 
ships for the properties of molecular clouds (Larson] |1981 j 



In particular, the mass-size relation reflects the hierarchical 
spatial structure of clouds ( jBeaumont et al. |2012 i. Molecu- 



lar clouds can agglomerate in small groups and chains, thus 
reflecting processes of the cloud formation: disintegration 
into smaller clumps or combining into bigger clouds. Due to 
the origin of the clouds from the large scale perturbations in 
the galactic disk molecular clouds can trace the large-scale 



structures, e.g. spiral arms and the bar (e.g. Wada & Koda 
2004). Several observational studies display existence of the 
galactic molecular ring (|Stecker et al. |1975|[Cohen fc Thad 



deus|1977 Roman-Duval et al. |2010 1 , however it is unclear 
whether this is an actual ring or just superposition of emis- 
sion from molecular clouds belonging to different spiral arms 
( Dobbs fc Burkert||20l2| ). 

In this paper we develop a model aiming to reproduce 
general characteristics of the formation of molecular clouds 
in the Galaxy, their physical properties, their distribution in 
the Milky Way using the 3D simulations with the compre- 
hensive gravitational potential, self-gravity of gaseous com- 
ponent, molecular hydrogen chemical kinetics, cooling and 
heating processes. 

The paper is organized as follows. Section 2 describes 
basic premises and equations of our model. Section 3 
presents numerical results on the disk evolution. Section 4 
presents analysis of the physical properties of the clouds in 
our model. Section 5 considers large-scale structures in the 
simulated Galaxy. Section 6 summarizes the results. 



2 MODEL 

2.1 Gas dynamics and galaxy potential 

The dynamics of the chemically reacting gas mixture in 
our model of the Galaxy can be written in the single-fluid 
approximation as follows: 

dp 



+ V ■ (pu) = 0, 



dt 

dpXi 
dt 

^ + V- (u® u) = -— - V(* Mt + * 9 ), 

dE 
~dt 



+ V ■ (puxi) = ps i: z = l,...,n s , 



Vp 



(1) 
(2) 
(3) 



+ V • ([E + p]u) = -u ■ V(* ea; t + 9 g ) - (A - T) , (4) 



where p is the gas density, p is the gas pressure, u = (u, v, w) 
is the gas velocity vector, n s is the number of components 
in the gas mixture (in our model n s = 2), Xi = Pi/p is the 

mass fraction of i-th component of gas, ^ Xi — 1j s i is the 

i=l 

formation/destruction rate of i-th component, E — p(e + 

^2 

— ) is the total energy, e is the specific internal energy, A 
and T are the cooling and heating rates, correspondingly, $ g 
is the gravitational potential of gas, ^ ex t is the total external 
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Figure 1. The rotation curve of the gas components in our model 
of the Galaxy, Vtot (thick solid line). The contributions from the 
stellar bulge, VJ,, the stellar disk, Vd, and the dark halo, Vh(r), ar e 
shown by dash, long dash and thin solid lines, correspondingly. 



gravitational potential produced by the massive dark halo 
^haio, the stellar bulge ^ bulge and the stellar disk ^disk- The 
latter includes the potentials of the spiral structure and the 
bar. The potential of the gaseous component is determined 
by the Poisson equation: 



A* 9 = 47rGp. 



(5) 



Hereafter we assume G = 1. The external gravitational po- 
tential can be written as follows: 



*ea;t = * halo + * 



bulge 



(6) 



Figure [T] present the rotation curve of the gas compo- 
nents adopted in our model of the Galaxy. For the maximum 
rotation velocity of the Milky Way we adopt a value ob- 
tained from the analysis of the recent observational data on 
the Galactic masers ( |Bobylev fc Bajkova |20l"o| ). The exter- 
nal gravitational potential of the Galactic subsystems, e.g. 
stellar bulge, stellar disk and dark halo, & sx t, produces this 
rotation curve with the parameters for the Galactic sub- 



systems fitted by Khoperskov & Tiurina (20031. The de- 



tailed description of the gravitational potential components 
is given in Appendix. 

In our calculations we use the value of the angular ve- 
locity Q. p = 31 km s _1 pc -2 as a fiducial value. This value is 
realized at the radius of corotation R c ~ 8 kpc (Bobylev & 
Bajkova |2010 1 . Our simulations have shown that variation 



of f2 p in the range 18 — 28 km s pc ( Amores et al. 



2009 



Sofue et al. |2009 Lepine et al. |2001 1 does not produce sig- 



nificant changes in the general properties of evolution of gas 
in the simulated Galaxy. 

2.2 Chemical kinetics 

H2 molecules in the interstellar medium are formed on 
the surface of the dust grains and dissociated by ultravi- 
olet Lyman- Werner photons and cosmic rays. We suppose 
that the dust density is proportional to the gas density be- 
cause the dust is well mixed with gas on the scales of our 
simulations. Thus, evolution of molecular hydrogen number 
density can be found from ( jBergin et fflT~||2004| ) 

dn2 
dt 



= ii 9r (T)nin 



[Ccr + (diss(N 2 , A v )]ri2, 



(7) 
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where n\ , ri2 are the number densities of atomic and molec- 
ular hydrogen, correspondingly, n — ni + 2n% is the total 
number density of hydrogen species, N2 is the H2 column 
density, R gr {T) = 2.2 x 10~ 18 ST 5 is the H 2 formation rate 
on dust grains (Tielens & Hollenbach|[T985), S = 0.3 is the 



efficiency of the H2 formation on dust ( Cazaux & Tielens 



20041, C 



6 x 10" 



s is the cosmic ray ionization rate, 



Ay is the extinction. Following Draine & Bertoldi ( 1996 1 



the photodissociation rate can be estimated as: 
C«K*.(#(HaMv) = C^ss{0)f S Meid{N{R 2 ))fdust, (8) 

where Cdiss(0) = 4.17 x 10 -11 sec -1 is the unshielded pho- 
todissociation rate, fshieid(N (H2)) is the H2 self-shielding 
factor, fdust is the dust absorption factor. 

To calculate the factors f shield and fdust we need to 
know both molecular, iV(H2), and total, N to t = iV(HI) + 
2iV(H2), hydrogen column densities. In order to do this ac- 
curately we should find a cumulative ultraviolet radiation 
field produced by the young stars and their clusters. Ob- 
viously, this problem is extremely complex, moreover, star 
formation processes are not included in our model. Hence, 



we follow a simplified approach introduced by Dobbs ( 2008 1 



We use a simple estimate that the column densities of the 
chemical species are just the local densities of these species 
multiplied by the typical distance to a young star, l p h'. 

(HI, H2) = lph ■ n(HI, H2). This distance is assumed to 
be a constant length scale for the whole disk. In our simu- 
lations we take l p h = 30 pc, which is in agreement with the 
number of massive stars expected from the Salpeter initial 
mass function. 

We assume that the Galactic gas has solar metallicity 
with the abundances given in |Asplund aT| ( [2005l ): [C/H] = 



2.45 x HT 4 , [O/H] = 4.57 x 10" 4 , [Si/H] = 3.24 x KT 5 . We 
assume that the dust depletion factors are equal to 0.72, 046 
and 0.2 for C, O and Si, correspondingly. Chemical kinetics 
of the heavy elements is not solved in our model the. We 
suppose that the carbon and silicon are singly ionized and 
oxygen stays neutral. This assumption is acceptable for the 
interstellar medium of the Milky Way due to existence of the 
strong ultraviolet radiation in the range ~10-13 eV (Habing 



19681 



2.3 Cooling and heating processes 

The cooling rates are computed separately for the tem- 
peratures below and above 2 x 10* K. In the low temperature 
range the cooling rate A in the energy equation |4| includes 
the typical processes of radiative losses for the interstellar 
medium: cooling due to recombination and collisional exci- 
tation and free-free emission of hydrogen ( Cen|199"2 |, molec- 
ular hydrogen cooling (Galli & Palla 1998), cooling in the 
fine structure and metastable transitions of carbon, oxygen 



and silicon ( Hollenbach fc McKee||l989 1, energy transfer in 
collisions with the dust particles ( |Wolfire et al. ||2003 \ and 
recombination cooling on the dust ( |Bakes fc Tielens||1994[ ). 
In the high temperature range, T > 2 x 10 K, the cooling 
rate for solar metallicity is taken from [Sutherlan d fc Dopita| 
( |1993[ ). We note that in our calculations the gas tempera- 
ture is generally below 10 4 K, higher temperature is reached 
only in small rarefied regions at the periphery of the disk. 

The heating rate F in the equation Q takes into ac- 
count photoelectric heating on the dust particles (Bakes & 



Tielens|1994]|Wolfire et al. |2003[ ), heating due t o H 2 forma- 
tion on the dust, and the H2 photodissociation (Hollenbach 
|fc McKee||1979 1 and the ionization heating by cosmic rays 
( |Goldsmitn et al. | 1978 1. 



2.4 Numerical methods 

To solve the system of gas dynamical equations ([l]|4| we 
use nonlinear finite- volume numerical scheme TVD MUSCL 
in the Cartesian coordinates (van Leer||1979| . The continu- 
ity equation for chemical species is constructed using the 
operator splitting scheme. Note that because in our model 
we have two chemical species, HI and H2, we can solve the 
kinetics equation for one of them only, e.g. for H2. At first we 
solve the continuity equation without right-hand terms and 
update the H2 density due to advection. Than we update 
this result after solution of chemical reactions. The equa- 
tion for H2 |7| is integrated using a fourth-order Runge- 
Kutta method. To find the gravitational potential of the gas 
|5]) we use the TreeCode method (Barnes & Hut 19861. 

The numerical resolution is 4096 x 4096 x 20 cells. The 
cells are assumed to be cubic with the sides of 7.3 pc. 
Therefore, the physical size of the computational grid is 
30 x 30 x 0. 15 kpc. We also made calculations with higher res- 
olution in the vertical direction, but lower in the disk plane 
direction, e.g. 1024 x 1024 x 100, and did not find significant 
difference in global characteristics of gas evolution. We ex- 
pect that this is sufficient to model the large scale dynamics 
of the Galactic disk and to study the formation of the spurs 
and molecular clouds in the vicinity of the spiral arms. 



3 DISK EVOLUTION 

Figure [2] shows the surface gas density at 50, 100, 200 
and 300 Myr. At the beginning the spiral shocks are formed 
due to the supersonic gas flow through the spiral gravita- 
tional potential of the stellar disk. After t — 50 Myr one 
finds a well-developed spiral structure in the gas disk (up- 
per left panel). The width of the gas spirals is significantly 
smaller than that of the stellar density wave due to higher 
stellar velocity dispersion. Owing to shear and wiggle insta- 
bilities the shock front becomes perturbed, and as a result 
the spurs and fragments are formed. Such spurs and frag- 
ments can be found at t = 50 Myr (upper left panel) and 
are clearly seen at t = 100 Myr (upper right panel). On one 
hand, higher density in such fragments stimulates molecular 
hydrogen formation and, as a consequence, efficient cool- 
ing of gas in the fragments. On the other hand, the cooling 
can lead to further fragmentation. We find that the efficient 
molecule formation starts at t ~ 30 — 50 Myr, which corre- 
sponds to the time scale of molecular hydrogen formation on 
the dust grains. A major part of the H2 molecules is formed 
in such fragments, which are the progenitors of molecular 
clouds. Note that a number of fragments formed may have 
super- Jeans masses. So, further fragmentation is expected 
to be amplified by both thermal instability and self-gravity. 

A large number of clouds is formed as a result of fur- 
ther evolution: a flocculent structure in the spiral arms can 
be clearly seen at t = 200 Myr (lower left panel of Figure [2j|. 
In the course of collisions the clouds can merge and form big- 
ger ones or can be disrupted completely. In the latter case 
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Figure 2. The evolution of the surface density of gas (g cm 2 ) at t = 50 Myr (upper left), t = 100 Myr (upper right), t = 200 Myr 
(lower left) and t = 300 Myr (lower right). The cross shows position of the Sun (8 kpc). 
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they produce a number of smaller cloudlets. The clouds mov- 
ing supersonically through the interstellar medium can be 
stripped due to the Kelvin-Helmholtz instability. Though H2 
molecules are formed and destroyed due to these processes, 
the total mass of molecular gas in the Galaxy increases until 
t ~ 200 Myr. Then it saturates because the quasi-station- 
ary regime is established between the molecule formation 
on dust and their destruction by the UV radiation from OB 
stars (lower left panel of Figure [2|. So, at t ~ 200 Myr one 
finds a population of clouds with sizes varying in the range 
~ 10 — 100 pc and typical lifetime about 10 Myr (in the 
next section we analyze the properties of the clouds in more 
detail) . 

At t = 200 Myr groups of clouds are settled along the 
spiral arms (Figure[2|. During further evolution these groups 
become larger both in size and mass. Gradually they stand 
apart from each other in the spiral arms. Such agglomer- 
ations of clouds have sizes more than 100 pc and the dis- 
tance between such associations at t — 300 Myr reaches 
several hundred parsecs. The agglomerations consist of sev- 
eral clouds with the typical number density £ 1 cm -3 . The 
typical thermal pressure p t h inside such clouds is about 
2-4 Mq pc -2 K which corresponds to the temperatures 
T ~ 80 — 200 K. Actual values of the gas temperature are 
lower because considerable portion of the internal energy is 
spent on the turbulent motions. However numerical resolu- 
tion of our simulations is not high enough to resolve the 
inner structure of the molecular clouds and the turbulent 
flows. The turbulence (especially small scale one) represents 
a common and complex problem. At this stage we are not 
able to extract the turbulent part of the internal energy. 



The velocity dispersion in a gas around molecular clouds 
is about 3-7 km s -1 and the velocity field around clouds 
shows a complex turbulent structure. To consider this we 
pick out a cloud with coordinates (-4.9 kpc, -4.9 kpc) at 
t — 300 Myr. Figure [3] shows velocity field around such 
cloud. The velocities are shown in the comoving coordi- 
nate system connected to the most dense part of the cloud. 
One can see the irregular form of the cloud and the collid- 
ing gaseous flows in the regions with enhanced density. The 
structure and physical parameters of the cloud locations ob- 
tained in our numerical simulations are close to those ob- 
served in the giant molecular clouds and their surroundings 
( Roman-Duval et al. 1 2010 1. 

Thermo-chemical processes taken into account in our 
simulations lead to formation of the two phases of the 
simulated gas with stable points around ~ 100 K and 
~ 10 4 K (Figure [ZJ see note on the value of the temperature 
above)). A gas with T ~ 100 K corresponds to the dense 
molecular clouds, i.e they have both high number density 
(n ~ 100 — 300 cm -3 ) and high H2 molecule abundance 
(x(H.2) ~ 0.3 — 0.5). These clouds concentrate to the spiral 
arms. A gas belonging to the warmer phase with T ~ 10 4 K 
is a rarefied atomic gas which resides mainly between the 
arms. A gas of the coldest phase (T ~ 10 K) occupies neg- 
ligibly small volume of space and is not considered here. 




x, [kpc] 

Figure 3. The surface gas density (g cm -2 ) and the velocity 
field (in the co- moving coordinate system) of the "cloud" with 
coordinates (-4.9 kpc, -4.9 kpc) at t = 300 Myr. 




Log(G) 

Figure 4. The phase diagram of thermal pressure p t ^ (MqK 
pc -2 ) versus surface density a (Mq pc -2 ) at the time t = 
300 Myr. The lines correspond to the constant temperature (see 
note on the value of the temperature in the text): 10 K (dashed 
line), 10 2 K (solid line) and 10 4 K (dash-dotted line). 



4 THE PHYSICAL PROPERTIES OF 
MOLECULAR CLOUDS 

In this section we consider physical and statistical prop- 
erties of the molecular gas in our simulations. This gas re- 
sides in the clumps which are called molecular clouds. A 
border of the cloud can be defined as location where H2 
abundance or H2 surface density exceeds given limit. The 
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Figure 5. Distribution of masses of the clouds for the H2 surface 
density thresholds S^ 2 = 1.5 X 10~ 5 g cm" 
10~ 4 g cm" 2 (dash line) at t = 300 Myr. 
to the empirical dependence N ~ Af~ 1,64 obtained by |Roman-| 
|Duval et al~\ \20W\ 



2 (solid line) and 1.5 X 
The line corresponds 



clouds in our analysis are delimited using a surface density 
threshold: over the 2D computational grid of the Galactic 
plane we find isolated groups of cells with the column den- 
sity exceeding given threshold. Such clumps normally have 
irregular form, and their effective linear size is estimated as 
a square root of their surface A. 

For our analysis we have chosen two values of the 
molecular hydrogen surface density threshold: Ef 2 = 1.5 x 
10~ 5 g cm -2 and 1.5 x 10 -4 g cm -2 . In both cases average 
density inside a cloud is greater than several dozens parti- 
cles per cc, which is close to the typical value of H2 number 
density obtained from observations of molecular clouds in 
the Galaxy (e.g. |Shu et al. [1987 1. 

We are using two thresholds for H2 surface density at 
the border of the cloud in order to show that our basic con- 
clusions on the statistical properties of clouds are indepen- 
dent of exact value of the threshold. This is important be- 
cause there is no single value of H2 surface density which 
exactly reproduces definition of a molecular cloud in obser- 
vations. 

Figure [5] presents the distribution of masses of the 



g cm 2 (dash line) 



both 

1.64. 



10~ b g cm (solid line) and 1.5 x 10 
at t = 300 Myr. 

We find that the spectra of cloud masses for 
thresholds follow the power law dependence N ~ Al -1 "': m 
the mass range M ~ 10 5 — 10 6 Mq for the higher threshold 
value and in mass range M ~ 10 6 - 10 7 M Q for the lower 
density threshold. This invariance of the spectrum on sur- 
face density limit indicates to that physical processes lead- 
ing to the formation of the structure with density threshold 
Ef 2 ~ 10~ 5 — 10~ 4 g cm -2 have the same nature. Hence, 
other statistical properties of the structures should be close. 
Note that the dependence N(M) for the higher threshold is 
valid in the same mass range as that obtained in the observa- 



tions of the Galactic molecular clouds (Roman-Duval et al. 



2010). For both thresholds one can see a lack of low-mass 



clouds. In the case of the higher limit the lack of clouds with 
M < 10 5 M© can be explained by our numerical resolution, 



o 10 6 - 




Size of cloud, [pc] 



Figure 6. Mass-size dependence for the H2 surface density 
thresholds Sf 2 = 1.5 X 10" 5 g cm (crosses) and 1.5 X 
10 -4 g cm~ 2 (circles) at t = 300 Myr. The lines correspond 
to several theoretical and empirical fits (M/Mq) = "/(A/pc 2 )^ 
given in |Beaumont et al. |l|2012|): (7, ft) = 240,0 .95 | |Larson|198"T) l, 
(7, ft) = 300, 1.5 jElmegreen &: Falgarone |1996] >, ( 7, ft) = 150,1.3 
<|Roman-Zuniga et al. \20W\ and (7, ft) = 228, 1*36 ( |Roman-Duval| 
\et al. |2010[ | shown by dotted, dashed, solid and dash-dotted lines, 
correspondingly. 



whereas similar deficiency in the observations comes from 
the incompleteness of the data ( |Roman-Duval et al. |20io| ). 
For the lower threshold the lack of clouds with M ^ 10° Mq 
can be explained by clustering of smaller and denser clumps. 
In other words, choosing lower surface density level we pick 
out more extended clusters of clumps, which contain both 
smaller and denser clumps. In the case of the lower thresh- 
old we cover larger regions with lower density, so that the 
increase of density limit should lead to decrease of size of 
the cluster and excludes extended low-density regions from 
consideration of. 

Figure [6] presents the dependence of mass of the cloud 
on its size (mass-size relationship) for the H2 surface den- 
sity thresholds Ef 2 = 1.5 x 10~ 5 g cm" 2 (crosses) and 
1.5 x 10" 4 g cm" 2 (circles) at t = 300 Myr. The sizes of 
the clouds vary between ~ 10 — 60 pc for the higher den- 
sity threshold. The most massive clouds have masses about 
2x 10 6 M Q . Thes e values are in good agreement with the ob- 
servational data ( Roman-Duval et al. | 2lnO). In some cases 
the sizes of molecular clouds can be as large as ~ 100—200 pc 
( jKirsanova et al. |2008[ ) , but these regions can be considered 
as a tight group of smaller clouds. In our analysis such ex- 
tended clouds can be found under the lower density thresh- 
old (crosses in Figure [6|: there are groups or chains of clouds 
seen in the periphery of the Galaxy at 300 Myr (Figure [2j|. 

The mass-size (or mass-area) relationship among molec- 
ular clouds, M ~ R 2 , which is well-known as Larson's third 
law, appears due to cloud evolution and numerous phys- 
ical processes in the interstellar medium. Many efforts to 
obtain or explain this universal relation are given in |Beau-] 
|mont et al. |p012| ), who tried to interpret the relationship 
in terms of the column density distribution function. How- 
ever, the scatter of the empirical fits is large (see lines in 
Figure [6] note that we pick out only several fits from the 
Table 1 presented in Beaumont et al. 20121. One can see 
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that the properties of clouds obtained in the numerical simu- 
lations are close to the power-law relation between the radius 
and mass of the cloud M = 228.R 2 ' 36 , which was obtained 
by applying a x-square minimization for the observed radii 



and masses of the Galactic molecular clouds ( Roman-Duval 



et al. 20101. There is a significant dispersion of our sim- 
ulated data around this fit. We find that our best fit of 
the simulated data M ~ E?' 4_2 - 16 (the ranges correspond 
different density thresholds) is close to that obtained in ob- 
servations. Note that our fit has a power-law index close to 
2 for both density thresholds considered here. 

Finally, we note that the velocity dispersion for the ma- 
jority of clouds varies in the range v t ~ 1 — 2 km s" 1 in our 
numerical simulations, that is a good agreement with the 
observational data ( Roman-Duval et al. pOlO l. 

Thus, we have found that the physical properties of clo- 
uds obtained from our simulation are close to those observed 
in the Galaxy. Our conclusions are coinstrained by the nu- 
merical resolution, because we can consider clouds with sizes 
larger than 10 pc, whereas even the parsec-size clouds are ob- 
served (Roman-Duval et aLj2010). Anyhow, our numerical 



results provide possibility to consider the global properties 
of molecular and atomic gas in the Galaxy. 



5 THE L — V DIAGRAM 

5.1 Molecular hydrogen 

Large-scale distribution of molecular gas in the Galaxy 
can be represented using the 'longitude-velocity' diagram 
(henceforth I — v diagram) constricted from observations 
of the CO(l-O) line emission ( |Dame et aF|[200T| ). Th e CO 



molecule is usually used as a tracer of molecular gas (Mai 



oney & Black 19881 because emission of much more abun- 



dant H2 molecule cannot be observed in the cold and warm 
regions. Current model of molecular kinetics in our simula- 
tions does not provide us with possibility to calculate abun- 
dance of CO molecule and the distribution of molecular gas 
is judged from the distribution of H2. Certainly, the CO- 
M2 conversion factor depends on many physical parameters 
( jShetty et al. |2011| [Feldmann et al. |2012[ ). However, there 
are no doubts that the general conclusion on agreement be- 
tween synthetic and observational distributions of molecular 
gas in the Galaxy can be drawn from comparison of I — v di- 
agram for H2 molecule in the model with CO data of [Dame] 
et al. I (12001 1. 



Figure |7| presents the I — v diagram for molecular gas 
at t = 300 Myr. The diagram is constructed for an observer 
located at the Sun position (see Figure [2j. One can find 
the large-scale structures corresponding to Perseus, Outer 
and Carina spiral arms. In the central part of the dia- 
gram, I ± (30 — 40)°, an intense and extended structure 
with the strong velocity gradient that passes through the 
origin is clearly seen. This structure by its location and gen- 
eral appearance resembles so-called Galactic Molecular Ring 



(Stecker et al. 1975 Cohen & Thaddeus 1977 Roman 



|Duval et al. ||2010[ ). However, Figure |2| shows that there 
is no pronounced ring of molecular material in our simu- 
lations. In our model the structure in the Molecular Ring 
locus of points of the I — v diagram is associated with the 
Galactic bar and the spiral arms at ~3-4 kpc. Thus, our 
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Figure 7. The l-v diagram for the molecular gas at t = 300 Myr. 
The mass of gas with a given (I, v) is indicated by grey scale: from 
small (light grey) to large (black) mass. 
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Figure 8. The l-v diagram for the atomic gas at t = 300 Myr. 
The mass of gas with a given (I, v) is indicates by grey scale: from 
small (light grey) to large (black) mass. 



results are in agreement with conclusions on the absence of 
the Molecular Ring drawn from the simple fitting technique 
Dobbs fc Burkert|(|2012|> and previous hydrodynamical sim- 



ulations (lEiiglmaicr & Gerhard 1999, Rodriguez-Fernandez 



fc Combe s|2008| |Baba et al. |2010| ). Note that the previous 



simulations did not include molecular chemical kinetics. 

We have to note on the existence of numerous small- 
scale structures in the diagram. They are obviously associ- 
ated with spurs, agglomerations of clouds and individual 
clouds. But we make no comment on that here because 
our consideration is concentrated on the extended structures 
scales of which greatly exceed our numerical resolution. 

5.2 Atomic hydrogen 

The atomic gas in our model is the main part of the 
gaseous component of the spiral arms and its distribution 
traces large-scale structures in the Galaxy. Figure [8] presents 
the synthetic I — v diagram for atomic gas at t — 300 Myr. 
The diagram is constructed for an observer located at the 
Sun position (see Figure[5|. The structures corresponding to 
the spiral arms are clearly seen in this Figure. The synthetic 
diagram does not show any prominent structure at the lo- 
cus of points corresponding to the Galactic Molecular Ring 
Figure |8j. Results of our modelling show that the atomic 
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gas in the simulated Galaxy is distributed much more uni- 
formly than the molecular gas. The synthetic I — v diagram 
shows reasonably good agreement with the observational 
data (see, e.g. I — v diagrams in |Kalberla & Dedes||2008| 



McClure-Grifiiths et al. 20041. So, our model reproduces 



main features of the observed distribution of the neutral gas 
in the Galaxy. 



"Dynasty" (Dmitry Zimin) for financial support. E.O.V. 
thanks to Russian Foundation for Basic Research (grant 11- 
02-01332). AMS was partly supported by the Russian federal 
task program Research and operations on priority directions 
of development of the science and technology complex of 
Russia for 2007-2012 (state contract 16.518.11.7074). The 
numerical simulations were made on the supercomputer fa- 
cilities of NIVC MSU "Lomonosov" and "Chebyshev". 



6 CONCLUSION 

In this paper we have studied formation of gas clouds 
in the model of our Galaxy. The 3D simulations have taken 
into account molecular hydrogen chemical kinetics, cool- 
ing and heating processes. Comprehensive gravitational po- 
tential have been used. It included contributions of self- 
gravitating gaseous component, stellar bulge, two and four 
armed spiral structure, stellar disk and dark halo. We have 
analyzed general properties of the simulated clouds and have 
compared them with statistical distributions taken from the 



recent surveys of molecular clouds (Roman-Duval et al 



2010) 



Our results can be summarized as follows: 

• the following stages of evolutionary sequence of the 
molecular cloud system formation were distinguished in our 
modelling: a) spurs and fragments (progenitors of molecu- 
lar clouds) are formed in the spiral arms due to shear and 
wiggle instabilities at the times t £ 50 Myr; b) numerous 
molecular clouds start to form after t XL 50 Myr which cor- 
responds to the time scale of molecular hydrogen formation 
on dust grains; c) a well-developed hierarchical structure of 
the clouds is formed at t ~ 200 Myr; d) at the moment 
of time t ~ 300 Myr the clouds agglomerate into extended 
associations with sizes exceeding 100 pc while the sizes of 
individual clouds vary in the range ~ 10 — 100 pc; 

• the total mass of molecules in the Galaxy increases 
until t ~ 200 Myr, when it saturates; the major part of 
molecule hydrogen is locked in molecular clouds where the 
H2 molecule abundance reaches x(H.2) ~ 0.3 — 0.5; 

• the statistical dependencies, e.g. the mass spectrum of 
clouds, the "mass-size" relation and the velocity dispersion 
of clouds, obtained in our simulations are close to those ob- 
served in the Galaxy ( Roman-Duval et al. |2010 l; 

• the structure is clearly seen in the Galactic Molecular 
Ring's locus of points on the simulated I — v diagram of 
molecular gas; this structure doesn't correspond to a real 
ring of molecular gas and likely arises due to superposition 
of emission from the galactic bar and the spiral arms at ~3- 



4 kpc, which supports conclusion made by ( Dobbs & Burk 
ert|2012 1 ; no prominent structure in the Galactic Molecular 
Ring locus of points is distinguished in the simulated I — v 
diagram of the HI gas. 
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APPENDIX A: EXTERNAL GRAVITATIONAL 
POTENTIAL MODEL 

The dark matter halo gravitational potential is taken in 

|199l| : 



the quasi-isothermal form (Begeman el 
M h 



C h 



ln(c) + ^ + l lni +f 



(Al) 



where £ = r/a^, Mh is the mass of the dark halo within a 
radius r^, ah is the scale of the halo and Ch = an^hjah — 
arctg(rh/ah)). In our model we take Th = 12 kpc, ah 
b kpc, M h = 0.64- 10 11 M Q . 

For the stellar bulge potential we adopt the King's 

r? ax (|Fridman & Khop- 



crskov 20111 



bulge 



rC b 



In 



r 

— + 



1+1^ 



(A2) 



with the following parameters Mb = 0.092 
0.2 kpc, rf M = 12 kpc, where 



10 n M Q , r b 



C b = ln(rr x /r b + y/l + (r™* /r b ) 2 ) 



(rr x /n)/(^l + (rr*/n) 2 ). 



For the three-dimensional stellar disk we take the ex- 
ponential distribution of surface density, then the stellar 



disk potential can be written as follows ( Binney & Tremaine 
|2008l 



^disk = TT<T z d ■ In (cosh (z/z d )) - 
-Tva r d -y-(I iy) K x (y) - h (y) K {y)) , 



(A3) 



where y = — - — , ctq = ^ d 2 i s ^ ne central surface density, 



2-?V 



27T-7-2 



M d = 0.4 • lO n M0, the radial scale of disk r d = 3 kpc, 
the vertical scale z d = 100 pc, Io(y),Ko(y),h(y),Ki(y) are 
the cylindrical Bessel functions of the first and second kind, 
respectively. 

Taking account the gravitational potentials of bar and 
spiral structure of stellar component we can obtain more re- 
alistic distribution of gas in the Galaxy. Following Cox & 
Gomez (20021 the perturbed potential of stellar disk 9^ iak 



can be written in the form of superposition of potentials gen- 
erated by the bar and the two and four-armed logarithmic 
spiral patterns: 



*d«fc = Vdish ■ (1 + £(*) ■ ^arms), 



(A4) 



where e(t) describes the evolution of the relative amplitude 
of the spiral stellar density wave. We assume a linear increase 

of the amplitude e(t < r) = 0.1 ■ - during r = 220 Myr, 

T 

after that time the amplitude is constant: e(t > r) = 0.1. 



10 Khoperskov et al. 

The relative perturbation ^arms is adopted in the form 
( Wada fc Koda|2004[ ): 



_ KlC0s(6l) K 2 COs(9 2 ) K 4 COs(6 4 ) , . . 

* arms — in i rTTn ' .1 /o 7 [A5 ) 

(1 + K1 ) 3/2 (1 + K2 ) 3/2 (1 + K4 ) 3/2 

where ki = {r/r a i) 2 , k 2 = {r/r a2 ) 2 , k 4 = (r/r ai ) 2 ; r al = 
2 kpc, r a 2 = 7 kpc, r a4 = 7 kpc are the radial scales for bar, 
two-armed and four-armed spiral patterns, respectively. The 
functions <di are: 

0i =2- 



6 2 = 2 •(</>- f 2 - cot(i 2 ) ■ In (r/r 02 )) , 

6 4 = 4 • (cj> - / 4 - cot(i 4 ) ■ In (r/r 04 )) , 

/i = / 2 = f& = 50° are the phases of bar and spiral patterns, 
ro2 = ro4 — 3.5 kpc, i 2 = i 4 = 18° are the scales and the 
pitch angles of two-armed and four-armed components. 



